Semi-relativistic effects in spin-1/2 quantum plasmas 
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Emerging possibilities for creating and studying novel plasma regimes, e.g. relativistic plasmas and dense 
systems, in a controlled laboratory environment also requires new modeling tools for such systems. This brings 
motivation for theoretical studies of the kinetic theory governing the dynamics of plasmas for which both rel- 
ativistic and quantum effects occur simultaneously. Here, we investigate relativistic corrections to the Pauli 
Hamiltonian in the context of a scalar kinetic theory for spin-1/2 quantum plasmas. In particular, we formulate 
a quantum kinetic theory that takes such effects as spin-orbit coupling and Zitterbewegung into account for the 
collective motion of electrons. We discuss the implications and possible applications of our findings. 
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I. INTRODUCTION 



Plasmas, in their full generality, make up a highly com- 
plex class of physical systems, from classical tenuous plas- 
mas, in e.g. fluorescent lighting, to dense, strongly coupled 
systems, such as QCD plasmas. The large span of plasma sys- 
tems implies that a wide variety of theoretical methods have 
been developed for their treatment. Even so, there are general 
principles that remain as common features between the differ- 
ent plasma systems. Therefore, methods used for one plasma 
system can in some cases be transferred to another plasma 
type, sometimes leading to new insights. One such exam- 
ple is the transferal of techniques for treating nonlinearities 
in classical plasmas to quantum mechanical plasmas. The lat- 
ter, often termed quantum plasmas (see, e.g., Refs. fl-IH]), to 
lowest order contains corrections due to the classical regime 
in terms of nonlocal terms, related to the tunneling aspects 
of the electron (in quantum plasmas, the ions are most often 
treated classically). Such tunneling effects can be incorpo- 
rated in both kinetic and fluid descriptions of the collective 
electron motion JUI^]. Such collective tunneling effects may, 
e.g., lead to nanoscale limitations in plasmonic devices |0] and 
bound states near moving test charges in plasmas JH]. Another 
mean-field effect that may be added to the dynamics of clas- 
sical plasmas concerns the electron spin, i.e. the possibility 
of large-scale magnetization of plasmas Jgt]. Such magneti- 
zation switching is known to be able to give new non-trivial 
features, such as metamaterial properties, allowing for, e.g., 
new soli ton modes flfl [Till . Moreover, the inclusion of the 
electron spin into the collective dynamics can be done either 
through a fluid or a kinetic approach fl2tl . Furthermore, the in- 
clusion of spin into the dynamics of a quantum plasma points 
in the direction of relativistic effects in such plasma systems, 
e.g. collective spin-orbit coupling. It is the intention of the 
present work to extend previous work into the weakly rela- 
tivistic regime. 

As indicated above, the dynamics of plasmas under extreme 
conditions is an important and integral part at many current 
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and up-coming experimental facilities, and such investigations 
therefore constitutes a highly active research field. In par- 
ticular, laboratory plasmas, such as laser generated plasmas, 
are currently presenting the possibility of studying previously 
unattainable plasma density regimes. It is well-known that in 
the high-density regime IU3I1 quantum effects start to play a 
role for, e.g., the dispersive properties of plasma waves 11141 — 
[lfill . In nature, such dense relativistic plasmas can be found 
in planetary interiors and in stars lfl7ll . Moreover, relativistic 
contributions to such plasma dynamics are under many cir- 
cumstances very important lfl8ll . Thus, when the parameters 
takes values characteristic for the quantum relativistic regime, 
one needs to consider more complex dynamical models in or- 
der to obtain accurate descriptions of a host of phenomena. A 
canonical starting point for dealing with high-density effects 
in a perturbative relativistic regime is offered by the quantum 
kinetic approach, here based on the Dirac description. Effects 
that can be included in such a perturbative model includes, 
e.g., spin dynamics, spin-orbit coupling, and Zitterbewegung. 
These examples have close connections to the nonperturba- 
tive relativistic quantum regime, in which e.g. pair production 
|[l9l - l22tl and other nonlinear quantum vacuum effects izl be- 
come pronounced 

Dense plasmas, and in particular short-time scale phenom- 
ena therein, have been successfully studied using Green's 
functions techniques, such as the Kadanoff-Baym kinetic 
equation s Izill (see also Refs. |[25l - l28tl for similar approaches, 
and Ref. B29I1 for an overview and examples from femtosecond 
laser physics). Although the Kadanoff-Baym equations and 
similar indeed gives ample opportunity to treat a wide variety 
of systems, their generality also makes simplifying assump- 
tions necessary, and under certain circumstances a mean-field 
model, that still retains memory effects and non-local struc- 
tures, can be an adequate approximation lfT7tl . In particular 
the mean-field approach is well suited outside the regime of 
strong coupling effects. Here, we will be interested in phe- 
nomena in plasmas that are not strongly coupled, but still in 
regimes where a classical plasma descriptions is not fully ad- 
equate. Here we stress that a large number of different dimen- 
sionless parameters (see e.g. Refs. lf30l [3lTl for more complete 
discussions) is needed to give a thorough description of var- 
ious plasma regimes. However, much insight can be gained 
by considering a simple density-temperature plot. In FigureQ] 
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Figure 1: Different plasma regimes in the temperature-density pa- 
rameter space. The dotted line is where the strong coupling param- 
eter r = Ep/ksT = 1, where E p is the potential energy due to the 
nearest neighbor. For higher densities the average kinetic energy 
of the particles is given by the Fermi energy kgTp rather than the 
thermal energy and the strong coupling parameter must there be re- 
placed by Tf = Ep/kg{T + Tf). In this graph Tf = 1 is illustrated 
by the dashed curve, and the strong coupling region (shaded) lies 
below this line. In this region higher order correlations has to be 
taken into account, and the mean-field description is not valid. For 
comparison, the lines TrCQp/ksT (dotted gray line), where hw p is the 
plasmon energy, and Tf /T (dotted-dashed gray line) are also plot- 
ted. These measure, respectively, the relative importance of wave 
function dispersion and the Fermi pressure. As a rough estimate, the 
quantum regime lies below these lines. The area marked ICF denote 
the regime of relevance for inertial confinement fusion experiments. 



a schematic view of the parameter re gime of interest for our 
study is presented, adopted from Ref. 13211 

In Section HH we recall the Foldy-Wouthuysen transforma- 
tion for particles in external fields and as a result we obtain 
a semi-relativistic Hamiltonian. We then go on to define the 
scalar quasi-distribution function for spin- 1/2 particles in Sec- 
tion [HI] In Section IV we then derive the corresponding evo- 
lution equation for the distribution function. The evolution 
equation clearly depicts the importance of the different terms 
of the relativistic expansion. A comparison of our results to 
previous studies is made, and we discuss the applicability of 
our equation as well as the interpretation of the variables in- 
volved. In Section V our theory is illustrated by means of two 
examples from linearized theory, and finally, in Section VI, 
our main conclusions are summarized. 



II. THE HIGH-ORDER CORRECTIONS 

The Dirac Hamiltonian can be written in the form 

H = Pmc 2 + Q. e + Q. , (1) 
where we have the even (fi£l e = £2 e j3) operator 

£l e = q§, (2) 



and the odd (j3£2„ = — i2 a j3) operator 

£l = cCC (p qA) 



(3) 



where j8 and Ct are the 4x4 Dirac matrices, m is the electron 
mass, q is the charge (q = — e for an electron), c is the speed of 
light, p is the momentum operator and <p and A are the scalar 
and vector potential, respectively. 

The odd operators in the Dirac Hamiltonian couples the 
positive and negative energy states of the Dirac bi-spinor. For 
the purpose of obtaining a perturbative expansion in the pa- 
rameter E/mc 2 , where E is the typical energy associated with 
the second and third term in (Q~|i, we assume that the first term 
in ([TJ is large compared to these terms. The consecutive ap- 
plication of the unitary Foldy-Wouthuysen (FW) lO trans- 
formation 



H -> exp (j5D. /2mc 2 ) Hexp (-j5Q. /2mc 2 



(4) 



yields a new Hamiltonian of the form (fTJ in which the new odd 
operators are of the order 1 jmc 2 . Performing this transforma- 
tion n times yields terms up to order (mc 2 )"". This gives a 
separation of the positive and negative energy states up to an 
arbitrary order n in \/mc 2 . 

Applying this transformation four times gives the following 
Hamiltonian for positive energy states with only even opera- 
tors H 



H = mc z - 
Hq 



-# + ^- fp--A 

Zm \ c 



4m 2 c 2 
1 

8m J c 



Ex 



qh 
2mc 
ih 2 q 



a B 



h 2 q 
8m 2 c 2 



a VxE 



V-E 



(5) 



where a here denotes a vector containing the 2x2 Pauli ma- 
trices, h is the reduced Planck constant, q is the charge, m is 
the mass, B and E are the magnetic and electric field and <p and 
A are the corresponding potentials in the Coulomb gauge. We 
see that the first four terms constitute the Pauli Hamiltonian, 
while the remaining terms are higher order corrections. In 
particular, the sixth and seventh terms together gives Thomas 
precession and spin-orbit coupling, while the fifth and eight 
terms are the Darwin term and the so called mass-velocity cor- 
rection term, respectively. 

In Eq. §5$ as well as in the Dirac theory we started from, 
the value of the spin g-factor is exactly 2. When applying the 
resulting theory, in section V, we will use the QED corrected 
value of g ~ 2.00232, however. In spite of the smallness of 
the modification it turns out that this correction is important, 
as the applications of our theory are very sensitive to the ex- 
act value of the g-factor. In fact, the sensitivity of the kinetic 
theory to the value of g was seen already in Ref. |35[| . This 
may suggest that for consistency, the Hamiltonian for QED- 
corrections should be added to Eq. ©. Such an approach 
would indeed modify the g-value to the desired one in the the- 
ory presented below, but the augmented Hamiltonian would 
also add several new terms in the evolution equation for the 
electrons. Those extra terms are at least smaller than those 
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kept by a factor of the order (g — 2), however. Thus the main 
effect from QED in the regime of study is the modification 
of the value of the ^-factor as compared to the Dirac theory. 
As a consequence, the contributions from QED (see e.g. Ref. 
l36ll for QED-corrections to the Dirac Hamiltonian) besides 
modifying the g-value will not be included here. 



III. GAUGE-INVARIANT STRATONOVICH-WIGNER 
FUNCTION 

The extended phase-space scalar kinetic model is obtained 
using the Hamiltonian ©. Following Ref. l32ll . we are able to 
construct a gauge invariant scalar kinetic theory using a den- 
sity matrix description for a spin- 1/2 particle. 

The basis states are |x, a) = |x) <E> \ where |x) is a state 
with position x and \a) is the state with spin-up a = 1 or 
spin-down a = 2. As a starting point of this model we use the 
spinor state y(x,a,t) = (x,a|i//) which fulfill the dynamical 
equation ihd t y/(x, a,t) = Hy/ (x,a,t), with the Hamiltonian 



With the spinors, we can define the density matrix as 
Pa/3 (x,y,0 = (x,«|p|y,j3) =£p ( ^ (x,a,/)^(y,j3,f) , 

(6) 

where p, is the probability to have a state i//,-. The density 
matrix fulfills the von Neumann equation 



ih 



d_p_ 

dt 



(7) 



Once the density matrix has been defined, we can define the 
Wigner-Stratonovich transform l38ll as 



W a p(x,p,t) = 



dh 
{2nti) 



■exp 



-z -4> 



Pa/3 



z z 

-;x— -,t 
2 2 



where the phase 



1/2 



4> = p--/ dr]A(x + T]z,t) 

c J-1/2 



(8) 



(9) 



is used to ensure gauge invariance of the resulting distribution 
function. The Wigner-Stratonovich transform has the property 
that it must to be taken separately for each component of the 
2-by-2 density matrix. 

Different approaches to construct a kinetic theory from 
the Wigner-Stratonovich transformation are discussed in Ref. 
J32tl . Following this reference we here define a scalar distri- 
bution function /(x,p,s,f) in the extended phase-space ll39ll 
where s is a vector of unit length. This distribution function 
satisfies that 

f(x,s,t) = J d 3 pf(x,p,s,t), (10) 

gives the probability to find the particle at position x with spin- 
up in the direction of s, and 

/(p,s,0 = J d 3 xf{x,p,s,t), (11) 



gives the probability to find the particle with momentum 
p with spin-up in the direction of s. Using the Wigner- 
Stratonovich transformation, the scalar distribution function 
will be defined as 0211 



1 2 

f(x,p,S,t) = — £ •<*)«/» lfy (X,p,t) 



a,p=\ 

— tr(l + s-<7)W(x,p,f), 

47Z 



(12) 



where tr denotes the trace over the spin indices. We recall that 
the expectation value polarization density is now given by 



(ff)(V)=tr[ffp(x,y,»)]=3 I d 3 pd 2 sf(x,p,s,t)s, (13) 



where we stress the need for the factor 3. This follows from 
the form of the transformation (fT2l and is needed to com- 
pensate for the quantum mechanical smearing of the distri- 
bution function in spin space. Furthermore, it should be 
stressed that the independent spin variable s constructed in 
(TTZl i generates the rest frame expression for the spin. In 
our theory, which is only weakly relativistic, this has lim- 
ited consequences. The relation between the rest frame spin 
s, and the spatial part of a the spin four-vector S is given 
by S = s+ [y 2 /(y+ l)](v-s)v/c 2 lE^l, where the kinematic 
quantities (i.e. the gamma factor 7 and the velocity v) can be 
expressed in terms of p and s (see below). Since our weakly 
relativistic theory presented here is only concerned with spin- 
dependent terms up to order v/c, the difference between S and 
s may be overlooked for the most part, e.g. when computing 
the magnetization current density. 



IV. EVOLUTION EQUATION FOR THE SCALAR 
DISTRIBUTION FUNCTION 



Using the above formalism we obtain a fully gauge invari- 
ant Vlasov-like evolution equation for charged particles. One 
of the most basic quantum effects is the tendency for the wave 
function to spread out. In the non-relativistic version of the 
theory [32j] this effect end up in operators sin(/iV r • V p ) and 
cos(KV x ■ Vp) acting on the fields and the distribution func- 
tion, where the operators can be defined through the trigono- 
metric Taylor-expansions l4lll . In our present theory we will 
view the spatial gradient operator V x as a small parameter, 
and drop terms of order or smaller, which means dropping 
particle dispersive effects, that are smaller by a factor of the 
order 8 2 where 8 is the characteristic de Broglie wavelength 
over the macroscopic scale length. The other approximation 
made is to only account for weakly relativistic effects, as de- 
scribed above. This implies that only terms up to first order in 
the velocity is kept, and that the gamma factor is put to unity. 
The evolution equation is found using the transformations ^ 
and ( fT2] i on the evolution equation (0, which together with 
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the above given approximations results in 



^ + {" + 2^ E ^ S + ^}- V ^ 



2ju / pxE\ 



pxE 



- ju (s + V f ) • <?M B 



2mc 



d'pf- 



h 2 q 
8m 2 c 2 



d<(V-E)d' p f, 
(14) 



where p is the momentum (which is related to the veloc- 
ity through the spin; see below) and ji = hq/2mc (or ji = 
ghq/Amc). 

The evolution equation ( fPfl ) has three new effects compared 
to the equation in Ref. IB2I1 for spin-1/2 particles. The first 
one is the Thomas precession effect where the previous the- 
ory 1I32I1 is extended by the substitution B->B-px E/2mc 
in the fourth and fifth terms of Eq. ( fPfl ) . This effect comes 
from the spin-orbit coupling contribution in Hamiltonian (0 
and, therefore, it is directly coupled with the evolution of the 
spin. The second new effect is the last term which is asso- 
ciated to the Darwin term. This term introduces the Zitter- 
bewegung effect of the electron, and is the only contribution 
proportional to h 2 . The third effect is seen in the velocity- 
momentum relation, which is highlighted in the second and 
third terms. In Eq. ( fPfl ) the term in {} brackets resembles a 
velocity which has been modified by the spin, which will be 
discussed in some detail below. Finally, we point out that the 
factor in front of V s f in the third term is indeed given by ds/dt 
l40ll . i.e. the laboratory rate of change of the rest frame value 
of the spin. Thus we note that Eq. (TBI) is consistent with the 
interpretation of s as the rest frame variable for the spin. 

Next we consider the relation between the velocity and mo- 
mentum. In order to relate this variables we use the Heisen- 
berg equation of motion for the velocity operator 



v = iM]. 

For the Hamiltonian (O we then get 



I(-/W-*A)-JUxE, 
m \ c J 2mc 



(15) 



(16) 



where we have neglected the last term in the Hamiltonian to 
simplify the equation slightly, (see the discussion about the 
mass correction below). We now recall that the Wigner trans- 
formation for an operator is multiplied by a factor (27z;/z) 3 , as 
compared to the the Wigner transformation for the density ma- 
trix. Similarly, for the spin transformation, the transformation 
for operators comes with a factor 3. Taking this into account 
and calculating the Wigner and Q transformation of the oper- 
ator above gives the final relation 



P 3u 

v = v(x,p,s,r) = - + - — E x s. 

m Lmc 



(17) 



This is the function in extended phase space, which can be 
used to calculate the average velocity and the current density 



of the plasma. An important question that arises is whether 
the current density based on the velocity (JT7J will give the free 
current density or if correspond to some other physical quan- 
tity, e.g. the total current density. This question is addressed 
below, where we calculate the energy conservation law of our 
system, which confirms that the velocity in Eq. (fTTI ) is indeed 
the variable corresponding to the free current density. 

In a more general context, the relationship between the mo- 
mentum and the velocity is nontrivial. E.g. the spin orbit cou- 
pling has been shown to arise as a Berry phase ter m 14211 . For 
further discussion of this interesting topic see e.g. II43I - I46I1 . 

When obtaining the evolution equation ( fT4b . we have not 
considered the effect of the mass-velocity correction term in 
order to get a more transparent formalism. This term will 
only produce a correction of the form p/m — > p/m(l + 
p 2 /2m 2 c 2 ) in the second term. Although this term is of the 
same order in an expansion in 1/c, as compared to other terms 
that have been kept, we will not consider it as the classical 
relativistic terms are already well-known. Instead we focus 
on the new effects introduced by the spin and the Zitterbewe- 
gung. 

The dynamics of the distribution function given by the 
Vlasov equation ( TT4b is in the mean field approximation cou- 
pled to the Maxwell equations in the form 

<9E 

V-E = 4np T , VxB=— +47Tj r , (18) 

at 

where the total charge density and total current density are 
given by 

<?P 

Pt-=Pf + V-P J r =J F +VxM + — . (19) 

at 

In the above expressions, the free charge density is 

p P = qJdSlf, (20) 

where d£l = d 3 vd 2 s is the integration measure performed over 
the three velocity variables and the two spin degrees of free- 
dom. The spin vector has a fixed unity length and it is thus 
convenient to use spherical coordinates (<p ( , 9 S ) do describe it. 
The free current density is given by 



[ 1 + 

m Lmc 



(21) 



With these charge and current densities, the conservation of 
charge is obtained as from ( TT4| ) to be d t pf + V • Jp — Q. Fur- 
thermore the magnetization M and the polarization P are both 
due to the spin and they are calculated respectively as 



M = 3jU / dtlsf 



and 



-3/f / dstp-f. 
lmc 



(22) 



(23) 
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The system of Maxwell's equations with the magnetization 
(l22l and polarization d23l and free current density (|2T1 i, to- 
gether with our main equation (TPfl) . satisfies an energy con- 
servation law of the form 

<?,W + V-K = 0. (24) 

Here the total energy density W is given by 

W = i(|E| 2 + |B| 2 )+|^(^|--3 J L i s-B^/, (25) 

and the energy flux vector K is given by 

K = Ex(B-M) + Jda(j^+3n ( R -^)' s ) v /- 

(26) 

Apparently the first terms in ( 1251 ) constitute the electromag- 
netic field energy density, and the integral term is the com- 
bined kinetic and magnetic dipole energy density. The first 
term of d26b is the Poynting vector, whereas the latter repre- 
sents the combined flux of kinetic energy density and mag- 
netic dipole energy. This energy conservation equation is a 



generalization of previous results for semi-classical theories 
for spin- 1/2 plasmas 113511 . It should be noted that although the 
theory presented here contains approximation, e.g. due to the 
weakly relativistic assumptions, the conservation law ( l24l is 
an exact property of the presented model. 



V. LINEARIZED THEORY 



In the present section, we are going to study the influence 
of the spin-orbit coupling and of the Darwin term on linear 
wave propagation. For this purpose we linearize the evolution 
equation dl4l . where the variables are separated into equilib- 
rium and perturbed quantities (using the subindices and 1, 
respectively, to denote them). Thus, the distribution function 
will be / = fo +fi, and the electric and magnetic field could 
be written as E = Ei and B = Bo + Bi respectively. The evo- 
lution equation to linear order becomes 



-£ + -• V -,/i + — P x B • V p fi + -£s x B • V,/i + nV xi [B • (s + V,)] V pi f, 
at m mc It 

= -qEi ■ V p fo - r^-Ej x (s + V,) • V x f — — p x Bi • V p f - p V,,- [Bj ■ (s + V,)] V pi fo 
2mc mc 

+ tt-V* [(P X Ei) • (s + V. s )] V pi fo - [Ei x (s + V. v )] X B • V„/ -^sxB r V s / 
Lmc Lmc L n 

+ X (P >< KO • VJb + ^V,(V • El )V,/ . (27) 



In the following, we only study electrostatic modes (e.g. 
Bi = in d27l i) propagating along Bo, as this gives a good il- 
lustration of the contribution from the relativistic terms, that 
are due to the Zitterbewegung effect and the spin-orbit cou- 
pling. 



A. Darwin term contribution 

Firstly we want to focus on the effect associated with Zitter- 
bewegung. The Zitterbewegung is a rapid oscillatory motion 
of the electron which implies that if an instantaneous mea- 
surement of its velocity is performed, the result is the speed of 
light. The amplitude of the oscillatory motion is x osc ~ h /2mc 
1134(1 , which means that the electron cannot be localized, but is 
rapidly oscillating in a volume of the order of the cube of the 
Compton wavelength. The Zitterbewegung is a quantum rela- 
tivistic effect and it is related to particle-antiparticle nature of 
the Dirac theory and to the nature of the spin. At the present, 
there is a growing interest in the detection of effects as the 



Zitterbewegung l37ll . 

The effect of the Zitterbewegung of the electron is intro- 
duced in the last term of Eq. d27l ), the Darwin contribution, 
which represents the smeared out electrostatic potential field 
that the electron sees when it fluctuates over a distance x osc . 

For the sake of simplicity, we examine the dispersion re- 
lation of Langmuir waves in an unmagnetized plasma. In 
order to focus on the effects of the Darwin contribution, we 
consider a one-dimensional unperturbed momentum distribu- 
tion and let fo — s- fo{p^)d(p x )d(p y ) (for a more realistic 3D- 
momentum distribution, even the electrostatic unmagnetized 
case couple to the spin terms, as we will see in the next sec- 
tion). The total distribution function will then have the form 
of f{z,p z ,t) = fa(p%) + fi(z,Pz,t)exp(ikz- ia>t), where f 
and f\ are the equilibrium and the perturbed distribution func- 
tions respectively. Furthermore, we consider a homogeneous 
plasma an neglect the motion of the heavy ions. The perturbed 
electric field is longitudinal, i.e. Ej = zE\ exp(/fe — /of). 
Using the evolution equation ( f27l ). the perturbed distribution 
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Figure 2: The parameter mvzitt/pF, relevant as a comparison when 
T — > 0, where pf is the Fermi momentum, is plotted as a function 
of density. This illustrates the relative importance of the Zitterbe- 
wegung and the Fermi momentum, i.e. the Fermi pressure, in the 
dispersion relation j3lT l above. 



function is then related to the electric field amplitude by 

-iqEi f_ h 2 k 2 



fi 



co — kp z /m 



I — — (28) 
8m 2 c 2 J dp z 



Combining Eq. (1281 1 with Poisson equation V • Ei = ikE\ = 
Anq f dQ.fi, we obtain the dispersion relation 



1 = 



k- 



h 2 k 2 



dp i 



./'() 



{p z Jm— co/k) 



(29) 



where the re-normalized distribution function /q 
fulfills/"^ dp z /o = 1. For phase velocities larger than 
than the characteristic spread in p z /m, we can Taylor expand 
the denominator, and write the dispersion relation as 



or = CO 




(30) 



where (p?) = f™„dp z p 2 fo is the average of the squared mo- 
mentum. Here the Landau damping term has been dropped, 
since the resonance is assumed to lie in the tail of the distribu- 
tion. 



The term proportional to h 2 k 2 is the Zitterbewegung contri- 
bution to electrostatic modes. For not too large wave-numbers 



it is more significant than quantum contributions from the 
Bohm potential [2] that scales as li 2 k 4 . Provided the wave- 
numbers are small, CO ~ C0 P , and the dispersion relation can be 
further approximated as 



C0~ 



(0 




(31) 



where vzitt = hco p /2mc can be understood as a velocity re- 
sponse of the plasma vzm ~ x osc co p to the rapid oscillations of 
the Zitterbewegung motion. The term v| itt A: 2 comes from the 
fact that the electron sees a smeared out electrostatic poten- 
tial, and therefore, a gradient of the electric field and a force 
of the order /i 2 V(V • E)/m 2 c 2 . Similar to the effects of the 
Fermi pressure (which gives a non-zero (/? 2 ) even when the 
temperature goes to zero) this implies a non-zero group veloc- 
ity of the electron plasma waves even in the cold case (see also 
Fig. 2 for a comparison between the Fermi statistics and the 
Zitterbewegung, relevant for Eq. (|3TT i when T — > 0). We note 
that the effect of the Darwin term related to Zitterbewegung 
becomes important for high-density plasmas when hco p /mc 2 
is not too much smaller than unity. 



B. Spin-orbit coupling contribution 

The spin-orbit effect of Hamiltonian Q appears as the 
Thomas precession correction (see e.g. Ref. 14711 ) of the mag- 
netic field in the fourth and fifth terms of Eq. ( IT41 ). 

Although this contribution can introduce interesting cor- 
rections to different types of wave modes, in this work we 
are going to follow the previous spirit and we analyze the 
quantum corrections to Langmuir waves, this time in a mag- 
netized plasma. These modes have been studied previously 
in Ref. ll48ll in a phenomenological relativistic formalism 
where the appropriate Thomas precession factor of 1 /2 was 
not used. We consider again longitudinal electrostatic modes 
Ei = £E\ exp(/fe — icot), but now propagating along an exter- 
nal magnetic field Bo = Bqz. As will be seen below, this will 
give an illustrative example of how the spin-orbit coupling 
modifies the usual dispersion relation. The distribution func- 
tion will be taken to be of the form /(z,p,s,f) = /o(p 2 , S ) + 
/i(z,p,s,?)exp(fe — icot). We use cylindrical coordinates for 
the momentum, i.e. p = xp± cos cp p + yp± sin cp± + zp z with 
p 2 = p 2 L + p 2 . Furthermore, we use spherical coordinates for 
the spin, i.e. s =isin0 v cosc/) s +j'sin0 5 sinc/) v +lcos0 s . 

At first order, taking the Fourier analysis of the evolution 
equation (f27), we have 
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p n d d \ r / /z 2 /t 2 \ df 



iknp ± Ex ( . n n d \ . . . .df 

sin Oi . + cos V -= — (cos (DpSinfflv — sin<p„cos <p s ) -=- — 

2mc \ d9 s J dp z 



U . . dfo 

- — p±E-[ (cos <p„ cos <p s + sin <p„ sin <p s ) 
Time o6 s 

qilBoEi ( . a a d \ ( df \ 

2mc2 I sm6 s cos <p s + cos 6 S cos (p s jQ- J \cos<p p j—\ 

^Imc^ ( Sin ° s Sin ^ + C ° S 01 Sin ^ ~dlT.) ( Sin ^ lfp~ ) ^ 



where fi) c = qBo/mc is the cyclotron frequency and co cg = 
(g/2)co c is the spin precession frequency [13311 - We note that 
the perturbed distribution function can be solved for in terms 
of the orthogonal eigenfunctions i//„ to the right hand side op- 
erator ll35l.l48ll . Accordingly we make the expansion 



1 



■ E E 8nrf(p±,Pz,9styn(p±,<Pp)exp(in'<p s ), 

\J2% „ = _oo„/ = _ 00 



where in general 

Vn{Px,(p P ) = 



1 



In 
1 



(33) 



exp[-/(n<p,, - (k ± p ± /co c m) sin (p^4) 



and /„« is the Bessel function. However, in this case for lon- 
gitudinal mode with k± = 0, then 



Wn(P±,<Pp) = —nF=e 



'<P P 



2% 



(36) 



Using the distribution function d33l l in Eq. d32l i. and then mul- 
tiplying both sides by y/*e~'"" f ' s and integrating over cp p and 
<p. v , we find that the only terms that survive in the sum d33l are 
gofi and g±i,=pi l32l [35l l48ll . Thus, we find the solution for f\ 



/1 



-iqE x {\ +h 2 k 2 /8m z c 



co — kp z /m 



ik\lp±E\ 
4mc 



■ '{<Pp-<Ps, 



2fimc 
e K<Ps-<Pp 



•K<Pp-9s. 



+ ■ 



co - Aco c - kp-J m co + Aco t - kp.J m 



co - Aco c - kp z /m co + Aco c - kp z J 



sin0 s + cos0. s — 1— ; 



dfo 
98, 



(37) 



where Aco c = co cg — co c . The expression d37b is combined with 

- icoEi = -\nJ Tz , (38) 

is used to deduce the dispersion relation. Due to the depen- 
dence on the angles cp s and cp p , the first term of (|37| | give raise 



to a free current density, whereas the other terms give raise to 
a polarization current density. The magnetization current den- 
sity vanish identically. Combining ( f37T > and d38l ) we find the 
dispersion relation 



1 



:+2 v zitt£ 



d£lp z 



-I — 

CO — kp z /m htn 2 c 2 



d£lp z 



A 



sin 0. 



1 



1 



sin0 4 +COS0J- 



dfo 



co - AcOc - kp z /m co + Aco c - kp z /m 
1 



d9 s J dp z \CO — Aco c ~kp z /m co + Aco c - kp z /m 



(39) 
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which is general for /o (where we have again used the distri- 
bution function re-normalized as / fodQ. = 1). 

As an example, let us examine an equilibrium distribution 
function with the form of a Maxwellian distribution and a spin 
dependent part Q2I1 



^'/^(l+cosft.) 
HT^^l-cos^) , (40) 



where T is the temperature and kg is the Boltz- 
mann constant, and the normalization factor is Nm = 
4n(7Zm 2 vf) 3 / 2 cosh(iiBo/k B T). We note that the expres- 
sion (|40b is the thermodynamic equilibrium distribution for a 



plasma of moderate density where the magnitude of the chem- 
ical potential is large l49ll . 



To simplify the integrals we will consider the frequency 
range where the wave frequency CO be close to resonance with 
Aco c . In this case, we neglect in (l39l l the terms with the 
denominators l/((3) + Afi) c — p z k/m) because they are small 
compared with the terms with denominators l/(o — Aco t — 
p z k/m). We are also going to take the limit when co — Aco t ^> 
p z kj m. Thus, using the equilibrium distribution function d40b . 
Taylor expanding the denominators, neglecting the poles in 
CO = p-k/m and co — Aco t = p z k/ m, and integrating over dD., 
we finally find the dispersion relation for Langmuir waves 
with spin-orbit coupling and Darwin effects 



or 



h 2 7t 2 co 2 p r 
8m 2 c 4 

1 



k 2 v} 



3k 4 vf 



1 Zitt 



(co-Aco c ) 2 
3k 2 v 2 



2(co-Aco c ) 4 



hn 2 co 2 vj 
Amc A 



tanhf^ 



k 2 v 2 



co-Aco c 2(co — Aco c y 



2co 2 



(41) 



The coefficient in front of the terms with denominators 
CO — Aco c are typically small, except for very strong mag- 
netic fields. Thus excluding the regime of an extremely 
strong external field, the frequency of the spin modes will 
be close to resonance, i.e. fulfill co ~ Ao c . More specif- 
ically, the deviation from exact resonance is of the order 
(co - Aco c ) / Aco c ~ (TiAco c /mc 2 )ta.nh(p.B Q /k B T). We note 
that spin induced modes with CO ~ Aco c have already been 
found by Ref. 113511 without the inclusion of spin-orbit cou- 
pling. However, it should be noted that the present wave mode 
is quite different from that previously found. In particular, the 
field is now completely electrostatic (whereas it was found 
to be completely electromagnetic in the previous case), and 
the present wave mode exists in the long wavelength regime, 
whereas the former wave mode ll35ll was dependent on a short 
wavelength, i.e. of the order of the electron gyroradius or 
shorter. Finally we note that even in the absence of an external 
magnetic field, a finite contribution from the electron spin re- 
mains, together with the Darwin contribution. In the absence 
of resonances we note that both these quantum contributions 
require a very high plasma density to be significant. 

VI. CONCLUSIONS 

The starting point for our theory is the FW-transformation, 
applied on the Dirac equation, which is used to pick out 
the positive energy states. Using the Wigner-Stratonovich 
transformation (HJ) on the density matrix, together with a Q- 
transformation M32I1 a scalar distribution function ( fT2b can be 
defined. In the weakly relativistic limit, the evolution equation 
is given by ( fT~4-b to leading order in the expansion parameters 



HL I p and jiB/mc 2 , where L is a characteristic macroscopic 
scale length, p a characteristic momentum of particles, and B 
a characteristic magnitude of the magnetic field. In addition 
to the magnetic dipole force and spin precession l32tl . which 
is included already in the Pauli Hamiltonian, Eq. ( TBI ) also 
contains spin-orbit interaction, including the Thomas factor, 
and also the contribution from the Darwin term. A compli- 
cation in the (weakly) relativistic theory that should be noted 
is the relation between momentum and velocity ( fT7] i, which 
now is dependent on the spin variable. Moreover, when clos- 
ing the system, the polarization ( l23l associated with the spin 
must be included in the current and charge density ( fT9b , and 
the expression for the free current density d2Tb is affected due 
to the aforementioned relation ( fTTI ). It should also be stressed 
that the spin variable used to define /(x,p,s,f) refers to the 
rest-frame spin, which is convenient, as two spherical angles 
for the spin variables is sufficient. 

In order to illustrate the theory, we have presented examples 
of electrostatic interaction in magnetized and non-magnetized 
plasmas. By picking a 1-D unperturbed distribution function, 
the influence of the Darwin term, which is associated with Zit- 
terbewegung, is highlighted in the case of a non-magnetized 
plasma. For a magnetized plasma, the spin-orbit terms leads to 
new types of resonances for electrostatic waves, which involve 
the combined effect of orbital and spin-precession motion. 

The theory presented here is of most interest for systems 
where at least one of the parameters hco p /mc 2 or jiB/mc 2 are 
not too small. Examples include e.g. laser-plasma interac- 
tions, astrophysical objects (e.g. white dwarf stars), solid state 
plasmas and strongly magnetized systems. The present paper 
is a first step to reach a fully relativistic quantum theory of 
plasmas. 
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